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SUMMARY 


A  multidimensional  descriptive  study  of  feline  lymphoma  was 
performed  using  age,  sex,  cell  type,  and  tumor  location  of  1733  cases. 

The  interactions  between  these  4  variables  were  determined  and  interpreted 
statistical!’.'.  Particular  consiaeration  was  given  to  the  relationships 
or  ace  to  tumor  location,  age  to  sex,  and  tumor  location  to  cell  type, 
which  were  shown  co  be  che  major  interactions  in  the  study.  The  use  of 
computer  analvsis  of  the  multiway  frequency  table  is  discussed.  The 
cuagestion  is  offered  chat  investigators  consider  the  use  of  the 
-.uitiwav  rrecuency  caoie  wnere  oossible  m  future  studies  of  various 
ciseases  to  ennance  the  depth  of  their  conclusions. 


Characteristics  within  a  disease  population  can  sometimes  help  to 
define  the  disease  itseif.  These  characteristics  or  variables  rarely, 
if  ever,  exist  without  some  mutual  association  to  each  other.  Employing 
these  associations  or  interactions  as  a  foundation  tor  a  biological 
hypothesis,  the  total  picture  becomes  more  complete  if  one  examines  all 
interactions  simultaneously.  Incorporating  many  variables  in  a  descriptive 
study  tends  to  emulate  nature  more  closely,  but  the  issue  of  explaining 
the  interactions  becomes  more  difficult  if  too  many  variables  are  used. 
Since  uata  in  general  have  a  quality  of  randomness,  it  follows  that  the 
mare  variables  examined  in  a  nopulation  study,  the  larger  the  sample 
sice  snouia  be  to  assure  adequate  representation  of  those  variables. 

Malignant  ivmphoma  is  considered  the  most  frequently  occurring 
cancer  in  cats,‘,+  ana  is  associated  with  a  virus,  feline  leukemia  virus 


The  age-related  incidence  rate  of  the  disease  is  bimodal,  peaking 

a 

at  around  1  and  10  years."  Males  have  a  2-fold  excess  in  incidence  over 

females.  Males  peak  in  incidence  at  about  7  years  of  age.  Females  peak 

at  aoout  i  year.  According  to  some,  malignant  Ivmphoma  has  certain 

3  4  10 

distinct  ceil  types.  ’  ,J"  Others  feel  that  the  various  types  are 
manifestations  of  the  same  cell,  but  simpiy  in  different  stages  of 
maturation.  J’  4 ’  *  From  the  standpoint  of  tumor  location,  the  internal 
form  tviscerai)  is  much  more  common  in  cats  than  the  multicentric  or 
solitary  forms.4 

A  method  of  simultaneous  examination  of  numerous  discrete  variables 

L  2 

exists  using  a  multiway  frequency  table.  ’  The  number  of  patients 
studied  must  be  at  least  10  times  the  number  of  cells  in  the  table. 

This  study  presents  such  an  ana Lysis  involving  age  groups,  ceil  type, 
sex  and  tumor  location  of  feline  malignant  lymphoma. 


wii 


Jateriaia  ana  Mecnoas  * 

The  stuuv  vas  conducted  using  the  records  of  che  Alameda/Contra 

Costa  Animal  Neoplasm  Registry.  A  complete  description  of  the  registry 

is  presentee  in  the  Liceracure.^  Tumor  submissions  were  chose  first 

-liagnosea  between  January  1970  ana  December  1978.  A  total  of  1728  cases 

vita  complete  information  was  submitted  during  that  period.  For  this 

jcuav,  ;ac:i  case  was  cacegorized  by  sex  (neutered  female,  intact  female, 

zeucerea  mate,  mcacc  male);  age  grout  (0-a  years,  5-3  years,  9  years  or 

4 

’ider,’;  rumor  location  (aoaominai.  tnoracic,  muiticencric ; ;  and  cell 

.  ~.r r.oovtia  veil— dir ferenciatea,  lynpnocytic  poorly-differentiated 

.'.istiocvric  .  .  The  age  grouos  represent  young,  mid-age  and  old  animals 

rusoeecivetv,  ;na  represent  the  two  peaxs  ana  che  trough  of  che  bimcdal 

o 

ige-rcisrca  mc-dence. '  The  following  coaes  were  used  to  designate  the 
var laoies :  A-age,  F-sex.  C-ceil  type,  and  L-cumor  iocacion. 

Cross— classification  of  the  1728  cases  ever  the  1  variables  created 
i  mux ci wav  ziaauency  taoie  (.Table  1).  It  consisted  of  108  cells,  5  of 
vniep.  a  a  a  zero  trequency,  prooaoiy  due  co  sampling  variations  and  would 
re  exoecteu  ro  ,;e  non-zero  if  the  samples  were  bigger.  [t  is  unlikely 
with  tn.e  icr:'i  sarauxe  size  in  this  study  that  any  empty  ceils  would  have 
a  true  expected  value  of  zero.  The  zero  frequencies  were  replaced  by  a 
frequency  of  one  co  facilitate  computer  analysis,  thereby  augmenting  the 
samoie  size  co  1733  cases. 

Analysis  of  the  muitiwav  frequency  was  done  using  a  program 

i  ”> 

designated.  3MDP3F."’“  This  particular  program  tests  for  associations 
among  cne  variables  (bivariate  and  more  compLex)  by  fitting  models 
involving  the  variables.  The  goai  of  this  type  of  analysis  is  to  find  a 
model  chat  has  minimal  complexity  while  exhibiting  good  fit,  which  is 


then  identified  as  che  "best"  model.  For  each  model  fit,  the  program 


2 

princs  a  goodness  of  tic  statistic,  G  ,  chd  log-likelihood  ratio.  This 
statistic,  for  large  samples,  is  equivalent  to  the  Pearson  chi-square 


statistic . 


calculated  as: 


all 


cens 


observed  x  log 
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observed 

expected 
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A  significantly  large  G  statistic  indicates  a  significant  discrepancy 
between  the  observed  frequencies  and  those  expected  for  that  model.  The 
program  has  the  capaoility  of  calculating  the  G”  value  for  all  models  of 
increasing  ccmoiexicv  starting  from  the  complete  independence  model, 

•  k)  <,$:  i.L)  in  tne  present  scuav.  to  the  completely  saturated  model, 

(ASuC ) .  As  tne  models  tested  become  more  complex,  the  G“  value  becomes 
smaller,  indicating  a  better  fit.  The  changes  in  G”  between  adjacent 
moaeLs  vita  corresponding  changes  in  degrees  of  freedom  are  used  to  test 
whether  increasing  the  complexity  of  the  model  achieves  a  significant 
improvement  in  the  goodness  of  fit. 

'.vhere  tne  uumoer  of  possible  models  is  large,  a  screening  procedure 
is  available  to  help  focus  on  the  several  models  that  might  be 
appropriate,  vniie  eliminating  those  modeis  that  warrant  no  consideration. 

-i 

The  program  automatically  princs  G“  statistics  for  the  complece 
independence  modeL  (A)(S)(L)(C)  in  the  present  study;  the  ail  possible 
first-order  interactions  model  (AS) (AL) (AC) (SL) (SC) (LC) ;  the  ail  possible 
second  order  interactions  model,  (ASL) (ASC) (ALC) (SLC) ;  and  so  forth  if 
higher  order  interactions  are  possible.  The  "best"  fitting  model  usually 
will  be  a  moaei  in  complexity  somewhere  between  the  most  complex  model 
automatically  fit  indicating  lack,  of  fit  and  the  least  complex  moaei 
automatically  fit  indicating  good  fit.  It  is  between  these  two  models 
that  subsequent  search  efforts  are  directed.  The  search  may  be 


accomplished  by  fitting  ail  models  in  complexity  between  these  two  models 
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?r  r-'  using  a  scenwise  selection  procedure.’  The  forward  stepwise 
pets  mure  starts  with  the  ieast  complex  model  automatically  fit 
indicating  lack  of  fit  ana  in  a  stepwise  manner  adds  higher  order 
interactions  until  a  good  fitting  model  is  achieved.  At  each  steD  the 
interaction  added  is  chat  which  produces  the  most  significant  reduction 
in  cr.e  sine  or  the  G“  obtained  for  the  modei  of  the  previous  step.  The 
selection  procedure  stops  at  that  steo  when  the  remaining  interactions 
:anr.oc  produce  a  statistically  significant  reauction  in  the  size  of  G~. 
The  oacxvara  elimination  orocecure  starts  with  the  most  complex  model 
i-csaatioai _v  tit  indicating  good  tit  and  in  a  stepwise  fashion  removes 
me  nignesc  order  interactions  until  a  model  is  fit  that  no  longer 
mo. laics  iooc.  tit.  At  eacn  scat*  che  interaction  removed  is  the 
mzeruccicn  mac  vnen  removed  makes  cne  least  impact  on  che  size  of  the 
resultant  i".  m  the  present  stud;.'  cne  screening  procedure  was  used 
wicn  :ocn  me  forward  stepwise  ana  backward  elimination  procedures. 
Jubsecuenciy  mi  possible  models  were  fit  as  a  check. 
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Sesuica 

The  initial  screening  showed  that  the  observed  data  were  not 
compatible  with  the  complete  independence  model  and  that  the  moael 
containing  ail  possible  first-order  interactions  gave  a  good  fit. 

Models  involving  second-  and  third-order  interactions  were  excluded  from 
further  consideration.  The  best  fitting  model  therefore  was  between 
(A)  (S)  CL)  (C)  ana  ( AS)  (AC)  (AL)(SL)  (SC)  (LG)  . 

A  slight  inconsistencv  emerged  in  the  forward  stepping  proceaure. 

The  searcn  went  through  5  stars,  yielding  two  models: 

•  AL)  (AS)  (.LC)  (SC)  (AC)  (1) 

and 

t’AL ;  (AS )  (LC )  ( AC )  (SL )  (2 ) 

These  two  meceis  were  both  indicated  as  suitable  because  the  change 
Ir.  .-t  stan  j  between  inclusion  of  (SC)  or  (SL)  was  very  close. 

Saint  the  backward  elimination  procedure,  the  selection  stopped 
after  (SC)  was  removed,  yielding  model  (2)  above.  In  light  of  this 
uncertainty,  it  would  be  justifiable  to  use  either  model  above,  or  the 
modal  with  ail  i  interactions: 

(AL )  (AS )  (LC )  (AC )  (SL )  (SC )  ( 3 ) 

1 

In  tne  forward  stepping  procedure,  major  reduction  in  G~  at  steo  1 
could  have  seen  made  by  including  anv  of  the  interactions,  (AL)  ,  (.AS)  , 
or  ’.LC)  in  tne  model.  In  fact,  they  were  the  first  3  interactions 
included  in  the  model  and  represent  the  three  strongest  interactions. 

The  remaining  first-order  interactions  were  of  a  lower  magnitude  as 
unown  by  tneir  lesser  impact  on  G~  at  step  L. 

Since  the  final  model  contained  nothing  higher  than  the  first  order 
(bivariate)  interactions,  further  statistical  interpretation  was  possible 
by  collapsing  the  i-way  table  (Table  1)  to  a  series  of  2-way  tables 
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rets  resent  mg  each  of  the  interactions  (Tables  2-7).  Each  cross 
ilassiz icatien  of  Tables  2-7  contains  2  numbers.  The  upper  number  is 
zr.e  ooservea  frequency  for  that  ceil.  The  lower  number  is  the  ratio  of 
;b served  frequency  to  the  expected  frequency. 

The  chi  sauare  values  for  Tables  2,  3,  and  4  were  above  90  and 
represent  cne  strongest  bivariate  raationshiDS  exhibited  in  the  data. 

,'^ble  1  describes  the  relationship  between  age  and  sex,  and  indicates 
tnat  vcung  intact  cats  of  both  sexes  occur  more  frequently  than  expected 
fr.  me  samoie.  Additionailv.  fewer  than  expected  intact  animals  of  both 
:exes  occurred  in  me  mid-ace  and  oid-aee  groups.  For  neutered  animals 
of  ootn  sexes  the  reverse  age  pattern  was  noted. 

Table  i  comoares  age  to  tumor  location.  In  young  cats,  there  was 
t  oaucitv  aodominai  cumors  and  an  abundance  of  thoracic  tumors.  In 
r. zc-  -na  old-age  cats,  tne  abdominai-thoracic  tumor  pattern  was  reversed. 
The  nuiticencric  tumors  showed  little  if  anv  age  predisposition. 

m  Table  1,  the  observed  number  of  muiticencric  and  abdominal 
tumors  with  tne  histiocytic  ceil  tvpe  was  higher  than  expected,  while 
the  veil-  me  pooriy-dirf erentiated  lymphocytic  ceil  types  in  multicentrii 
.-,nd  .lococmai  locations  were  observed  to  be  fewer  than  expected.  Well- 
uid  ooor l-’-u if f erentiated  lymphocytic  thoracic  tumors  were  abundant 
while  histiocy cic  tumors  of  the  thorax  were  sparse. 

Although  the  chi  square  values  for  Tables  2-7  were  all  significant, 
those  for  Tables  2,  3  and  4  were  considerably  higher  and  the 
observed  :exoected  ratios  for  Tables  2,  3  and  4  generally  had  a  wider 
variation  from  a  value  of  one.  Nearly  aLi  the  observed  frequencies  in 
Table  5  (cell  type  by  sex  group)  were  close  to  the  expected  frequencies. 
The  exceotion  to  this  was  the  paucity  of  histiocytic  cell  types  among 
intact  femaLes  and  of  lymphocytic  well-differentiated  3mong  neutered 
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males,  as  well  as  a  mild  abundance  of  lymphocytic  poorly-differentiated 
among  intact  females. 

Table  6  is  somewhat  difficult  to  summarize.  Suffice  it  to  say  that 
among  young  animals  there  was  a  paucity  of  the  histiocytic  cell  type  and 
that  among  the  mid-age  and  old-age  groups  the  histiocytic  tumors  were 
more  abundant  than  expected. 

In  Table  7,  it  appears  significant  chat  the  "female"  row  provided 
most  of  the  departures  from  expected  frequencies;  that  is,  more  thoracic 
and  less  abdominal  and  multicentric  tumors  than  expected. 

Discussion 

Some  of  the  reasons  for  conducting  a  study  such  as  this  are  to 
define  certain  associations  between  disease  characteristics  that  may 
suggest  interactions  of  various  biological  phenomenon;  or  to  point  out 
a  possible  similarity  with  results  of  studies  conducted  from  a  different 
prospective,  but  arriving  at  the  same  point.  Additionally,  this  studv 
illustrates  the  use  of  the  multiway  frequency  table  analysis,  a 
procedure  not  too  often  used  in  medical  research. 

A  second-order  interaction  is  one  in  which  the  pattern  of 
association  between  2  variables  is  dependent  on  the  level  of  a  third 
variable.  One  can  perceive  from  this  how  a  third-order  interaction  may 
behave.  It  is  unfortunate  that  only  first-order  interactions  emerged 
from  this  study,  but  the  method  of  analysis  is  illustrated.  Without 
this  method,  it  would  be  tempting  to  conclude  that  trivariate  and  higher 
interactions  do  exist,  based  on  only  visual  inspection  of  the  tables. 

Such  a  conclusion  would  be  done  without  the  benefit  of  any  statistical 
interpretation  such  as  was  attempted  here. 


%  4. 
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3cudicis  have  shown  that  age  plays  a  part  in  initiation  of  the  ■ 

13-16 

disease.  )c  Che  incacc  animals  in  our  study,  perhaps  natural 

selection  through  lack  of  resistance  eliminates  them  from  the  population 
sc  a  voung  age,  leaving  fewer  than  expected  in  the  older  groups.  On  the 
jt.;er  aana.  since  the  reverse  -  as  true  of  neutered  animals,  perhaps  sex 
.-.or-.ones  pis-.-  1  part  in  expression  of  the  disease,  once  infected  with 
cr.e  virus. 

lir.'.sisr icy  and  possibly  an  explanation  of  the  paucity  of  histiocytic 
vail  tvoes  cr.  tnoracic  tumors  might  be  found  in  a  study  of  human 
~.s__  mar.  t  .•v.pucma  patients  in  Saudi  Arabia*'  in  vnich  there  was  a 
reneenc  ::r  nistiocytic  ceil  types  to  present  as  abdominal  lesions. 

i  n 
Q 

Anocnar  :v.:man  studv*  indicates  chat  a  predominance  of  one  B  cell  line 
la  aSoOc:::iU  vicn  weil-uif ferentiated  ivmpnocvtic  ceil  type  and  cases 
vim  :  accro  isea  3  ceil  population  are  primarily  of  the  poorly- 
;if f ■‘rer.tiatcu  lympnocytic  ceil  tvpe.  This  same  report  suggests  that 
future  morpno logic  studies  of  the  disease  should  include  simultaneous 
m.mui'.Oi'gic  classification.  Using  the  multiway  frequency  table,  one 
couiu  oMD.mu  cnac  suggestion  ip  include  tumor  location,  which  was  shown 
to  interact  strongly  with  ceil  type  m  our  study.  Some  say  thoracic 

19_?1 

cumocs  ,,re  primarily  T  ceils  and  abdominal  tumors  are  3  cells. 

Although  this  studv  was  intended  mainiv  to  further  characterize  the 
disease,  additional  studies  using  the  same  techniques  but  with  different 
variables  may  shed  light  on  prognosis  and  therapy  of  the  disease. 
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relationship  of  age  to  sex  (AS) 


Sex 


Intact 

niale 

Neutered 

tale 

Intact 

temaie 

Neutered 

female 

132 

547 

207 

233 

!  * 

0.340 

1.233 

0.940 

■,2 

195 

32 

72 

3.763 

1.397 

0.567 

0.360 

40 

12  S 

32 

97 

0.675 

1 . 3U3 

0.675 

1 . 380 

ji  with  o  d. : . 

p  <  0 , OUi . 

an  frequency/ expected  frequency. 
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TA3LS  3 — 3ivariate  relations'niD  of  age  to  tumor  location  (AL) 


Location 

Young 

Age 

Mid-age 

Old 

Mult Ltencria 

Observed 

393 

136 

108 

Ratio 

1.000 

1.025 

.970 

Abdominal 

Observed 

353 

169 

172 

Ratio 

0.824 

1.169 

00 

■H 

rhoraic 

Observed 

323 

56 

23 

Ratio 

1.303 

0.669 

0.327 

Obi  acuare  = 

=  102.339  vich  -4  d.f. 

o  <  0.001. 

Ratio :  Observed  crequeucv/ expected  frequency. 
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TA3LZ  5 — Bivariate  relationship  of  sex  Co  cell  Cvpe  (SC) 


Ceil  Cvpe 

Sex 

Intact 

Male 

Neucered 

male 

Intact 

Female 

Neutered 

female 

LvT.pnoc vc  1c .  well- 

Observed 

103 

L50 

33 

111 

differentiated 

Racio 

1.024 

0. 363 

1.187 

1.070 

■iisciocvcia 

Observed 

177 

32S 

93 

136 

Ratio 

1.003 

1.082 

0.759 

1.023 

Lymonocvric.  poorly- 

Observed 

110 

192 

95 

105 

differentiated 

Ratio 

0.974 

0.989 

1.210 

0.902 

Chi  s auare  =  10.1205 

vich  6  d.f. 

p  v  0.005 

Ratio :  Observed  frequency/ expected  frequency. 


